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Abstract 

Food webs with intraguild predation (IGP) are widespread in natural habitats. Their adaptation and 
resilience behaviour is principal for understanding restructuring of ecological communities. In spite of the 
importance of IGP food webs their behaviour even for the simplest 3-species systems has not been fully 
explored. One fundamental question is how an increase of diversity of the lowest trophic level impacts the 
persistence of higher trophic levels in IGP relationships. We analyze a 3-species food web model with a 
heterogeneous resources and IGP. The model consists of two predators directly coupled via IGP relation and 
indirectly via competition for resource. The resource is subdivided into distinct subpopulations. Individuals 
in the subpopulations are grazed at different rates by the predators. We consider two models: an IGP 
module with immobilization by the top predator and an IGP module with species turnover. We examine the 
effect of increasing enrichment and varying immobilization (resource transfer) rate on a stable coexistence 
of predators and resources. We explore how the predictions from the basic 3-species model are altered 
when the IGP module is extended to multiple resource subpopulations. We investigate which parameters 
support a robust coexistence in the IGP system. For the case of multiple subpopulations of the resource 
we present a numerical comparison of the percentage of food webs with stable coexistence for different 
dimensionalities of the resource community. At low immobilization (transfer) rates our model predicts a 
stable 3-species coexistence only for intermediate enrichment meanwhile at high rates a large set of stable 
equilibrium configurations is found for high enrichment as well. 

Keywords: Intraguild predation, Immobilization, Alternative resource, Multiple resource traits, Stable 
coexistence 



1. Introduction 



In spite of the pr evalence and importance 



of omnivory food webs ([Pimm and Lawton . 



1978 



Vadeboncoeur et all 120051) in natural communities 
their population dynamics to date remain poorly 
understood, even for only three species in the 
community. Even in simple systems a plethora 
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of nonlinear effects such as flexible consumer be- 



haviour ( Leibold et al 



2005), intraspecific interac- 



tions betwe e n com peting; consumers and resources 



(Holt et al 



1994), inhomogeneity of the env iron- 



ment (j Amarasekare , 



adap tive foraging 



2007 



Janssen et al 



Krivan , 



1996; 



2007) and 



Krivan and Diehl 



20051 ) precludes easy theoretical treatment and in- 
terpretation. 

One example of a non-trivial omniv ory food web 



is a system with intraguild predation (IPohs et al 



1989 



Finke and Denno 



2002 



Borer et al 



20031) 



Intraguild predation assumes that the same or- 
ganism is both competitor and predator to an- 
other member of the food web. The IGP mod- 



eluding coexistence ( 
alternative stable sta 


Polis and Holt 




1992 


and 


tes ( 


Holt and Huxel 




2007 



Daugherty et al 



2000 



2007). Simple mathematical 



models (Polis and Holt 



Namba et al 



1992 



Diehl and Feifiel 



20081 ) have been evoked in at- 



tempt to explain the persistence of IGP interac- 
tions in natural habitats. However predictions from 
the mathematical theory of 3-species IGP systems 
state that a high resource carrying capacity pro- 
motes the exclusion of intermedi ate trophic levels 



and thus destabilizes interactions ( Diehl and Feifiel 



20011 ). What is puzzling that various empirical 
studies of omnivory document however coexistence, 
but not exclusion, over the entire range o f nat 



ural resour c e pro ductivities ([Mvlius et al 



Borer et al 



2001 



2003). On the basis of experimental 



observations a theore t ical 3- s pecies omnivory model 



I Stoeckcr and Evans, 



Diehl and Feifiel 



1985 



Holt and Polis 



1997 



for the communal resource (jDiehl and Feifiel 120011 ) . 
Yet empirical data suggest a rob ust persistence 



of IGP systems in both terrestr ial (jBrodeur et al 



2000: 



Arim and Marquet 



munities (Polis et al 



Borer et al 



2003; 



19891 



2004) and aquat ic com- 



Mvlius et al 



Denno and Fagan 



2003) 



2001 



Theoretical models that are focused on the 
aspects of stability and coexist ence of species 



in 3-level systems with the IGP ( Polis 



1992 



Holt and Polis 



1997 



Abrams et al 



nd Holt 



1994 



20101 ). as a rule, largely reduce the complex- 



i ty of interactions ob served in realistic systems 



(Thomson et al 



200711 . Such oversimplifications 
can influence the population dynamics as well as 
critically impact species persistence. Even though 



three species ( 


Polis and Holt . 


1992; 


Holt and Polis. 


1997 




Diehl and Feifiel 




2000 


2001 


) a number of 



empirical studies deal with larger food webs that 
involve more than three species potential l y en 



gaged in IGP interactio ns (jRosenheim et al 



Woodward et al 



20051) . 



1993 



Spatiotemporal heterogeneity of the environ- 
ment often is invoked as one of the explana- 
tory mechanisms for the coexistence between mul- 
tiple species com peting for the same resources 



( Hutchinson 



1961). It has been observed that such 



a spatiotemporal heterogenei ty can affect the d i 



2006) 



versity in prey populations ([Amarasekare , 
Indeed an inhomogeneity in prey items that share 
common resource and predators is critical in deter- 
mining the responses of ecological community. For 
systems with multiple prey composition various co- 



200 ll) predicts the coexistence existence patterns can be fou nd depending o n the 



only at superior competitive abilities of the IG prey levels of resource productivity (jLeibold 



1996). It is 



not clear yet how the diversity in a prey community 
will affect the behaviour in the IGP systems. 

The effect of a habitat structure on the IGP is 



discussed in various recent models ( Amarasekare, 



2006 



2007 



Janssen et al 



2007). For example a 



stable coexistence of the intraguild prey due to 
inhomogeneity of a habitat can be supported by 
creating temporal refuges for prey and reducing 
the encounter rates among preys and predators 



( Janssen et al 



20071) . In addition the stability of 



the IGP can be enhanced by an inclusion of addi- 
tional factors such as behaviourally mediated effects 



( Janssen et al 



20071) . 



To include the effect of an increasing diversity of 
resource and IG predators on popu lation dynamics 
recen tly the 3-species IGP model (jHolt and Polis . 
19971 ) was modified by Holt and Huxel (2007). 
The authors extended the basic 3-species omnivory 
model to the so called " partial IGP" model in which 
" partial" overlap among competitors for a single re- 
source exists and both predators have exclusive re- 



sourc es to exploit. It was shown (jHolt and Huxel 



20071 ) that an alternative resource enhances the 
tolerance of the IG prey against attacks from IG 
predators. Independently of a competitive status 
of the IG prey in exploitation for a shared re- 
source it can persists by utilizing an alternative re- 
source. An extended formulation of the IGP model 
with trophic supplementation has been proposed by 
Daugherty et al. (2007). The authors investigated 
three forms of a supplementary feeding outside of 
the basic IGP module and postulated a higher po- 
tential for persistence of the IG prey due to its ef- 
ficient exploitation of external resources. 

There is growing evidence that in many sys- 



tems the IG prey has a mutualistic or at least 
facilitative relationship with the IG predator 



(|Crowlev and Cox . 



20111 ). Including such facilita- 



tion in ecological theory will fundamentally change 
many basic predictions and will enable a better un- 
dersta nding of functioning of many natural commu- 



nities (Bruno et al 



2003). 



Especially in the IGP systems an emphasis 
should be given to the elucidation of the effects of 
faci litation on community c omposition and stabil- 
ity ([Crowley and Coxl . 120111 ) . Contrary to the com- 
petitive exclusion principle in systems with com- 
petitors for a s ingle resource stabil ity stems from 



commensalism (Hosac k et al 



2009). Hereby one 



consumer can in some way alter the habitat to ben- 
efit the other. Recently such an interaction was 
observed in experime nts with a microzo oplankton 



food web community (jLoder et al 



xxxxi ). The ex- 



perimental system included two predators: a tintin- 
nid species Favella ehrenbergii and a heterotrophic 
dinoflagellate species Gyrodinium dominans. They 
are both grazing on a phototrophic dinoflagellate 
Scrippsiella trochoidea. The authors showed that 
the IG predator F. ehrenbergii can precondition a 
substantial part of the common resource S. tro- 
choidea during its feeding procedure by immobiliz- 
ing the common prey without ingestion. Such pre- 
conditioned individuals can be captured more easily 
by the IG prey G. dominans than the mobile indi- 
viduals of the same resource species. This mutual- 
istic interaction leads to higher growth rates of the 
IG prey in the presence of the IG predator. The au- 
thors characterized their experimental observations 
as a facilitative IGP relationship with a commen- 
salistic pattern. Our motivation for this modeling 



study was to investigate if such commensalistic pat- 
terns can create loopholes for a stable coexistence 
of all species in the investigated system. Of our 
major interest was if in the IGP system an immo- 



bilization (Loderetal 



xxxxi ) or the partitioning 



of prey populations into distinct groups of individ- 
uals offers opportunities for competition avoidance 
among both consumer species. 



posed in ( 


Polis and 


lolt 




1992 




Holt and Polis 


1997 


Diehl and Feifiel, 


2000 


2001 


) to include mul- 



tiple subpopulations of prey. Furthermore, we ex- 
plored the effect of diversification of the resource 
available to higher level consumers on the species 
persistence by numerical simulations of an extended 
IGP module. Specifically, we investigate how the 
addition of new links to a focal IGP module en- 
hances stability of population dynamics by reducing 
the competitive interactions of predators for their 
shared resource. 

In order to explain the results of the experimen- 
tal findings of Locdcr ct al. we investigated the 
influence of multiple traits of the resource commu- 
nity on a stable coexistence in the 3-species model 
with different types of resource. For this purpose 
we adapted and reformulated the original model 
by Holt and Polis (1997) and added a new type of 
interaction. This link specifies the immobilization 
mechanism that depends on the densities of mobile 
and immobile resource items and the top predator 
which creates the immobile resource fraction during 
feeding. The immobilization term is used to model 
the interactions between the IG predator and the 
resource. 

Another type of interaction considered in this 



paper is a resource turnover mechanism. This 
mechanism describes mutual interactions between 
species from distinct resource subpopulations. The 
interaction term depends exclusively on the re- 
source subpopulation densities. The rate of 
turnover is constant. If no turnover or immobi- 
lization of individuals from one group to another 
occurs then the basic IGP model with a single pop- 
ulation of resource is recovered. We discuss the 
influence of immobilization and transfer of species 
on the coexistence patterns in a system with dif- 
ferent subpopulations of the resource and compare 
the results with the basic 3-species IGP. 

This paper is organized as follows: in the first 
section we introduce a general 3-species IGP model 
with a new type of interaction that links the re- 
source pools to the top consumers. In the following 
sections two distinct IGP formulations with n = 2 
resource subpopulations are discussed. Both mod- 
els are derived from the basic IGP module by in- 
cluding additional links: (i) the immobilization by 
the predator and (ii) the resource turnover. In the 
Results section we numerically investigate stability 
of equilibrium densities for various trophic configu- 
rations. Data from numerical analysis are presented 
for the IGP model with the immobilization and for 
the model with the resource turnover. At last we 
discuss results for a general IGP model with the re- 
source turnover mechanism and n > 2 subpopula- 
tions of the resource. After sketching the main con- 
clusions we review the model predictions and com- 
pare t heir relevance to t he immobilization experi- 



ment (Loderetal 



xxxxi ). Furthermore we discuss 



possible alternative reformulations of the model. In 
the Appendix explicit forms for the steady states for 



two simple analytical cases and multidimensional 
system are specified. As a part of a linear stabil- 
ity routine the Jacobian matrices for two types of 
formulations are given. Finally, we carry over to 
a higher dimensional formulation and describe the 
parameters choice and the equilibrium densities. 

2. General model 

We introduce an omnivory model with an IGP 
unit derived from a simple non-spatial Lotka- 
Volterra system with the linear functional responses 
adapted from Holt and Polis (1997). The origi- 
nal model consists of populations of two predators 
(IG predator and IG prey) and a common resource. 
Here, we include new features such as a resource 
differentiation mechanism which affects palatability 
of a fraction of resource for the predators. Specif- 
ically, the entire resource population is subdivided 
into distinct groups under the assumption that the 
groups differ from each other by the quality and fit- 
ness of the individuals. They are consumed by the 
predators at different group-specific grazing rates. 
The differentiation of the resource could be due to 
damage by the predator or initial inhomogeneous 
distribution of the resource quality. Afterwards, we 
generalize our model to the case of the multiple re- 
sources. 

The food web model for a multiple number 
of prey subpopulations {Sk}k—i lS sketched in 
Fig-Hall] The top predator F and the intermediate 
predator G are engaged in the IGP and share a com- 
mon resource Si. The resource pools are not inde- 
pendent because there is an exchange of individuals 



among different subpopulations {Sk}u—2 following 
the links in Fig.QJi. Another special case of the IGP 
with two distinct populations of resource S m and Sj 
is presented in Fig. [Tp. Shown is a schematic view 
of trophic interactions including intraguild preda- 
tion and two populations Si and S m of immobilized 
and mobile resources respectively. The IG prey G 
competes with the IG predator F for both resource 
types and is also an additional resource for the IG 
predator. The size of the population Si increases 
due to immobilization of individuals from the pop- 
ulation S m by the IG predator F. 



We begin with an overview of a general IGP 
model and all the important trophic links and pa- 
rameters that are used to define it. Later we focus 
specifically on two different formulations of the gen- 
eral IGP model. 



The general model for a food web with an in- 
homogeneous resource is d erived from the Lotka- 



2000 



Volte rra omnivory model (|Diehl and Feifiel , 
20011 ) with the interaction term that accounts for 
the transitions among different pools. The Lotka- 
Volterra omnivory model consists of n + 2 equa- 
tions. It is used as an approximation for the food 
web community with the IGP and n > 2 mutually 
interacting subpopulations of the resources. In the 
absence of predation a basal po pulation Si devel 



ops a ccording to logistic growth (jDiehl and Feifiel 



20001 ) . The set of equations for the population den- 
sities are written as follows: 



ifiere and everywhere in the text the numerical subscripts denote species at the same trophic level. 

5 



Shared 1st resource : 



dSi 



^± = [ r (\ - SxK' 1 ) - aG - fF]S x 
dt L v ; J 

— zi(S\, S2, • ■ • , S n , G, F), 



Shared kth resource (k — 2 . . . n) : 
dS k 



dt 



= Wk{Si,Sz...,S n ,G,F) 
-[b k G + f k F + m k }S k , 



Intermediate predator (IG prey): 



dG 
~dt 



Z2\Sl, S*2, • ■ • , S n , G, F) 

-(gF + m g )G, 



Top predator (IG predator): 



dF 
~dt 



Z3\S\, S-2, ■ ■ ■ 1 S n , G, F) 

+(g'gG-m f )F, (1) 



The parameters of the model and main popula- 
tions are described in details in Table [TJ Here r is 
the maximum specific growth rate of the resource 
population Si, K is the carrying capacity of the 
resource defin ed as enrichment facto r in th e pre- 



vious models ( Diehl and F eificl 



2000 



2001 ). The 



sub-populations {S k } k=2 are derived from the basal 
resource Si via immobilization or via individual-to- 
individual turnover. Species from Si and {S k }k^i 
are consumed by the IG predator at potentially dif- 
ferent rates / and {fk} k =2 and by the IG prey 
at rates a and {bk} k= 2 respectively. The differ- 
entiation among subpopulations {S k }k = 2 is P re ~ 
served by a choice of distinct predation pressures 
{M£=2,{Mfc=2, feeding rates {f k }% =2 , {b' k YU 
and mortality coefficients {wfc}£ =2 . The density- 
independent mortality rates for Si , G and F are 
mi,m g and m/ correspondingly. They are used as 



factors limiting the growth of the populations in 

A key assumption of the model is that there is 
only one-directional movement between the basal 
resource Si and its fractions {Sk}k^i- The local in- 
teractions among individuals from alternative pools 
are embedded via functional terms {wk}k=2 P r0 ~ 
vided in Table [2] for each type of the IGP formu- 
lation. These terms account for transitions among 
the resource items {Sk}k=ii- The general omnivory 
model (TTJ can be reduced to three types of IGP 
formulations: system with immobilization and sys- 
tems with the resource turnover for n — 2 sub- 
populations and for n > 2 pools. For each of 
the formulations specific expressions of functional 
forms zi, £2,^3 and {wk} are provided in the Ta- 
ble [2j The term zi is responsible for the exchange 
of individuals among subpopulations {Sk} due to 
the species turnover or the immobilization mecha- 
nism. The transfer of individuals from the popula- 
tion Si to {Sk}k=2 happens instantaneously at con- 
stant rates {ck}k =2 correspondingly. Analogously 
{qk,j}k^j are defined as instantaneous migration 
rates among subpopulations {Sk}k= 2 - The terms 
Z2 and Z3 are used to evaluate the total predation 
of the IG prey and the IG predator on the resource. 

To achieve a stable persistence of all species the 
IG prey should benefit more from an alternative 
resource than the IG predator. For this reason, 
whereas the attack rates of the IG predator are 
equal for different resource pools, the IG prey es- 
tablishes a higher predation pressure on subpopu- 
lations {Sk}k^i than on the basal pool Si. The 
numerical values for the attack rates are chosen 
to be close to the experimentally observed values 



(jLoder et all Ixxxxl ) . 



Holt and Huxel (2007) used an extended IGP 
module with alternative resources that are defined 
independently. They evolve according to their own 
intrinsic growth rates. As opposed to the formu- 
lation given by Holt and Huxel (2007) and to a 



mode l with trophic supplements (jDaughertv et al 



20071 ) here we do not consider external alternative 
resources. In our model with immobilization the 
population density in every resource pool varies due 
to immobilization by the IG predator and consump- 
tion by the predators. Similarly in the formulation 
with the resource turnover the transfer mechanism 
between resource subpopulations plays a role in ex- 
change among the distinct resource pools. Alter- 
native pools grow due to the influx of species from 
the basal resource or the other pools. Therefore 
the sizes of subpopulations are controlled mainly by 
a number of direct encounters with the IG preda- 
tor (immobilization) or by a species turnover from 
one resource subpopulation to another. In addition, 
the individuals in the different pools of the basal 
resource are distinguished by group-specific preda- 
tion pressures that establish a top-down regulation 



of densities of each subpopulation. 

In the following sections we present an explicit 
formulation of the model with immobilization and 
of the model with the resource turnover for n = 2 
subpopulations. 






s 2 








Figure 1: (a) General structure of the food web model with 
two predators and multiple resource subpopulations; (b) the 
structure of the food web with the presence of the immo- 
bilization mechanism by the IG predator. The resource is 
subdivided into populations of mobile S m and immobilized 
Si individuals. The links represent: (solid) food resources 
for the top predator, (dashed) food resources for the inter- 
mediate predator, (dot-dashed) transitions between different 
resource pools. 
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Table 1: The variables and parameters for the general model (fTJ. 



Definition 



General model 

Populations 



s k 


population size of resource in the £;th pool, 


G 


population size of IG prey 


F 


population size of IG predator 


Parameters 




r 


maximum specific growth rate of the resource S\ 


K 


resource carrying capacity or enrichment 


a 


attack rate of predator G on Si subpopulation 


f 


attack rate of predator F on 5*1 


fJ 


attack rate of predator F on G 


Ck 


per capita effect of species S\ on Sk 


Ik, j 


per capita effect of species Sk to Sj 


bk 


attack rate of predator G on Sk species 


fk 


attack rate of predator F on Sk species 


m k 


mortality rate of species from Sk subpopulation 


m g 


mortality rate of G 


m f 


mortality rate of F 


9' 


converting efficiency of food resource G into F 


Si 


growth rate of F from resource Sj 


b 'j 


growth rate of G from resource Sj 
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Table 2: The variables and parameters for the models J2J and lt3t . 



Definition 



Model with immobilization 

Populations 

S m , Si population sizes of mobile (immobilized) resource 
Parameters 

r maximum specific growth rate of population S m 

K resource carrying capacity 

a, b attack rates of predator G on mobile (immobilized) population 

/ attack rate of predator F on mobile and immobilized populations 

i m immobilization rate 

a' , /' conversion efficiency factors 

Model with a resource turnover 

Populations 

S\ , S2 population sizes of resources 
Parameters 

r maximum specific growth rate of subpopulation Si 

K resource carrying capacity 

a, b attack rates of predator G on subpopulations Si, S2 

f attack rate of predator F on subpopulations Si and S2 

t r transfer rate or per capita effect of Si on 5*2 
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Table 3: Description of the functional forms used in the system JTJ. 



Description 



Model equations 



System with immobilization: 



Si = S m , S 2 = Si,z\ = i m S m F, w 2 = zt, z 2 = a'(aS m + bSi)G, 

Z3 = f'f(S m + Si)F 



System with resource turnover: z\ = t r SiS 2 , w 2 — z\, z 2 = a'(aSi + bS 2 )G, z 3 = /'/(Si + S 2 )F 
(n = 2 resource subpopulations) 

System with the resource turnover: z\ = S\ Y^k= 2 c kSk, z 2 — G^™ =1 b'jSj,z 3 = fjSj, 

(n > 2 resource subpopulations) Wk = Sk{ckSi + Xw=2 Qk,jSj), k = 2 . . .n 



2.1. System with immobilization by predator 

The system with the immobilization illustrated 
in Fig.[T]b is derived from the equations ([T]) for two 
resource subpopulations by substituting the inter- 
action terms Z\,z 2 , Z3 and w 2 from Table [TJ After 
the substitution the set of equations for the IGP 
model with immobilization yields: 
Mobile resource: 

dS n 



dt 



= [r(l - S^nK- 1 ) - aG]S n 
-(/ + i m )FS m , 



Immobilized resource: 
dSi 



= i m FS m -{bG + fF]Si, 



IG prey: 



dG r 1 

— = [a ab n 
dt 



a'bSi - gF - m g ]G, 



IG predator: 
dF 



dt 



- [ffiSi + S^ + g'gG-m^F, (2) 



where the state variables S m and Si are the densi- 
ties of mobile and immobilized species. Note that 
the feeding rates of the top predator F on both 
populations Si and S m are equal. By contrast, the 
attack rate of the IG prey on immobilized subpopu- 
lation is higher than on mobile species. The relation 
b > a holds in the presence and in the absence of 
the predator F. This assumption is well justified 
by the observati ons of an experimen t with artificial 



immobilization ( Loderetal 



xxxxf ). G. dominans 



demonstrate a strongly selective behaviour towards 
immobilized species when offered in a mixture with 
mobile cells of S. trochoidea. It was measured that 
ingestion rates of the predator in the immobilized 
prey treatment were by a factor of 20 greater than 
those in the control treatment. 

The stability of equilibrium densities and the 
persistence zones of the system ([2]) with a non-zero 
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immobilization rate are discussed in Section 13.11 

2.2. System with the resource turnover 
2.2.1. General case of n = 2 subpopulations 

The model with the resource turnover is derived 
from the general case (JXJ) by substituting the func- 
tional forms from Table [5J It is written as follows: 
1st resource: 



dSi 



[r(l - SiK- 1 ) - oG]5i 

-[fF + trS^Sx, 



2nd resource: 



^ = [t r S 1 -bG-fF]S 2 , 



IG prey: 

dG 
~dt 

IG predator. 

dF 
~dt 



[a'aSi + a'bS2 — gF — m g ]G, 



[f'f(S 1 + S 2 )+g'gG-m f }F (3) 



All the parameters are chosen the same as for the 
system with immobilization (J2J . Note that the evo- 
lution equations are written as in © but immo- 
bilization term is replaced with the transfer term 
that is dependent on the population densities. The 
transfer between the two subpopulations occurs 
each time whenever species from two different pools 
encounter each other. In the simplest case the num- 
ber of encounters is proportional to the population 
densities of Si and S 2 . 

If the density of second subpopulation is zero 
and no differentiation in the resource takes place at 
t r = than the top predator F outcompetes the 
predator G due to a higher predation rate (/ > a). 
This outcome is predicted by the basic IGP model 



(jDiehl and Feiflell 12000 , 120011) . By contrast, when- 
ever the turnover of species takes place and non- 
zero densities are produced in the resource pool S 2 
the intraguild predation introduces a higher pres- 
sure on the second subpopulation S 2 . This will 
potentially lead to a negative effect on the popu- 
lation density in S 2 and to higher levels of subpop- 
ulation Si . The result of this interaction is that the 
3-species coexistence is reached via the IGP com- 
petition trade-off. 

3. Main results 

We illustrate an emergent dynamical behaviour 
for the three formulations provided in Table [3] with 
stability diagrams. Due to high dimensionality of 
the models CE}-© the analysis of an entire param- 
eter space is intractable. Only several illustrative 
examples for every formulation will be shown here. 

3.1. Model with immobilization 

In Fig. [2] the regions of stable positive equilib- 
rium solutions versus immobilization and enrich- 
ment are shown. The parameter space is parti- 
tioned into several stability zones associated with 
the regions of coexistence, exclusion of both preda- 
tors and exclusion of the IG prey at G = 0. The 
boundaries defined for partitioning of the diagram 
are found from the eigenvalue analysis (see Ap- 
pendix). As shown in Fig. [5] at low enrichment the 
densities of both predators decay to zero and the 
summed abundance of the resource reaches steady 
state at S m + Si = K. The case of zero im- 
mobilization Jiasjjej3njdread^_^ 



ous studies ( Holt and Polis . 



1997 



Diehl and Feifiel 
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2000). At low immobilization and at high enrich- 
ment only the top predator and resource are stable 
and positive, just as in the 3-species IGP model 



(jDiehl and Feifiel . 



20001 ) . whereas the coexistence 
between both predators and common resource is 
possible only in the regions of intermediate enrich- 
ment. A higher mortality rate for the predator G 
results in its extinction in the region of low immo- 
bilization in Fig. ^Bp where an extra resource can no 
longer support its persistence. Only the IG preda- 
tor and resource persist in this region of parame- 
ters. Situation is different for higher immobilization 
where a large region of coexistence for both preda- 
tors exists. The equilibrium densities shown on the 
diagrams are defined in the eq. (|A.5[) in Appendix. 

Fig. [3] shows equilibrium densities of the four 
components of the food web and their dependence 
on enrichment and immobilization rates. For a high 
immobilization rate the resource population is dom- 
inated by immobilized individuals. Meanwhile at 
low immobilization mobile and immobilized popu- 
lations increase along the gradient of enrichment 
an adverse pattern occurs at high immobilization. 
The growth rate of the IG predator is noticeably 
reduced at i m — 0.2 due to an increase of the com- 
petitive trade-off with the IG prey. The dependence 
of the population densities on the enrichment of re- 
source is shown in Fig. [4] for the immobilization 
i m = 0.3. Meanw hile as predicted from t he stan- 



dard IGP model (Diehl and Feifiel 



20001) the IG 



prey is excluded at high enrichment in the model 
with immobilization at i m = 0.3 the IG prey ben- 
efits from immobilized resource and its persistence 
is increased at a broader range of carrying capac- 
ities. The density of the mobile (immobilized) re- 



source subpopulation reach saturation threshold at 
a higher enrichment (see Fig. U]). 
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Figure 2: Regions of stable coexistence for the immobi- 
lization model J2}. Stability diagrams are partitioned into 
the regions of stable (unstable) coexistence and alterna- 
tive states with exclusion of one of the predators. Let- 
ters stand for persistence of different trophic configurations: 
(S) =exclusion of both predators, (Si+Sm+F) =coexistence 
of top predator and resource, (Si+Sm+F+G) =3-species co- 
existence. Parameters are: r = 0.4, / = 0.12, a = 0.025, b = 
0.1, a' = 0.8,g = 0.025,/' = 0.2, g' = 0.5, 771/ = 0.04. Two 



plots for: (a) 



0.02, (b) m g = 0.06 are shown. 
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Figure 3: Equilibrium biomass densities along immobiliza- 
tion and enrichment gradients. Shown are densities of: (a) 
mobile resource, (b) immobilized resource, (c) predator G,(d) 
predator F. Parameters are used as in Fig. [2] 
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Figure 4: Equilibrium biomass densities versus enrichment 
for fixed immobilization rate i m = 0.3. 



Figure 5: Regions of stable coexistence for the immobiliza- 
tion model J2J. The level curves for the grazing rate a de- 
fine the boundaries of the stability regions. The color bar 
shows the values of a. Parameters are: r = 0.5, / = 0.1, b = 
0.1, a' = 0.8,g = 0.07,/' = 0.2, g' = 0.5, m g = 0.02, m f = 
0.04 

How sensitive is a stable coexistence to small 
variations of the attack rates of the intermediate 
predator G? Will our predictions be still valid? To 
examine the system behaviour for different attack 
rates of G the coexistence zones are exemplified for 
different values of a in Fig.[SJ Colorcode is assigned 
according to grazing rate a. Overall the stability di- 
agram exhibits similar pattern as in Fig. [2] Specif- 
ically, the region of stable coexistence enlarges for 
higher immobilization. As it seems reasonable the 
number of stable solutions and the 3-species per- 
manence zone in Fig. [5] gradually broadens with the 
increase of predation pressure from predator G. Si- 
multaneously fewer exclusion steady states for the 
predator G are discovered. 
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3.2. Model with a resource turnover 

3.2.1. Case of n — 2 subpopulations 

In this section the equilibrium solutions and sta- 
bility of the equilibria are discussed for the sys- 
tem ((5]) with the mechanism of species turnover. 
In Fig. [6] the regions of stable (unstable) equilibria 
are plotted versus the enrichment and the transfer 
rate. The results are contrasted on the stability di- 
agrams in Figs. \6\ a — d for different predation rates 
of the IG predators. Four different states are lo- 
calized in the parameter space that corresponds to 
stable (unstable) persistence and exclusion of the 
IG prey (IG predator). We investigate how the dy- 
namics in the extended IGP system responded to 
variation of enrichment levels. For each case shown 
in Fig. [6] an increase of enrichment is accompanied 
with a series of bifurcations in the system mani- 
fested by an invasion of higher trophic levels simi- 
lar to predictions from the linear food chain theory 



( Oksansen et al 



19811 ). For instance, at low en- 



richment both regimes 4 and 1 are stable. Further 
increase of K at fixed t r = 0.05 results in a chain of 
bifurcations from a stable regime 4 to an unstable 
2 and subsequently to a stable coexistence regime 
1. A further increase of carrying capacity favours 
an exclusion of G and shifts the population densi- 
ties towards the dominance of the IG predator. An 
interesting feature is that at low transfer rate only 
a coexistence of the IG predator and the resource 
is found. The second subpopulation S2 is extin- 
guished fast due to predation and low transfer rate. 
The steady states found for low enrichment are sim- 
ilar to the case of a single prey population without 
transfer mechanism at t r — and S2 = 0. As typ- 



ified on the diagram in Fig. [5] d the IG prey levels 
remain positive. Since the IG prey has an advan- 
tage as a competitor for the shared resource only 
the IG predator gets excluded from the system. 

The stability behaviour of the system ^ is 
highly sensitive to the alternations of attack rates of 
G and the productivity of resource r. Changes of 
these parameters produce different emergent pat- 
terns as shown in Fig. [JJ The location of states 
of stable (unstable) permanence and the exclusion 
zone of G is still comparable to the patterns shown 
in Fig. |51 however the region of 3-species coexis- 
tence gets visibly reduced. The reduction is more 
evident on the plots Fig. [Jj a, c and d. At higher 
transfer rates the coexistence of all 3-species is no 
longer observed and only the population of IG prey 
and resource persist. Due to low productivity the 
densities of the basal resource Si are quickly de- 
pleted and the IG predator is driven to extinction. 
On the contrary, conditions become more profitable 
for the IG prey that is released from the IGP pres- 
sure and simultaneously obtains more benefits by 
predation on the extra resource S2. 

At low transfer rates (Fig. [7J; and d) the IG 
predator is excluded independently on carrying ca- 
pacity of the resource. As expected, with increase 
of the attack rate of G the population of the IG 
predator is driven to extinction due competition 
with IG prey. However, situation becomes more 
favourable for the IG predator at higher values of 
the transfer coefficient t r . For high enrichment and 
intermediate transfer the IG prey is excluded from 
the system. At a fixed enrichment several alter- 
nating states are found along the gradient of t r 
(see Fig. [Jj d). For example, at K > 2.5 the be- 



ll 



haviour of the food web is very sensitive even to a 
small alternations of t r . Indeed, the system passes 
through distinct steady states just within a small 
increment of transfer rate. The exclusion of the IG 
predator is observed at t r < 0.05, the coexistence 
is found at t r 0.06 and the exclusion of the IG prey 
is achieved at t r 0.07. Finally at a higher transfer 
values (t r > 0.14) both predators enter the system 
and persistence is reached. 

After presenting the results for the systems (0) 
and ([3]) we proceed to a more complex situation 
with n > 2 of distinct subpopulations of the re- 
source. 

3.2.2. Case ofn>2 prey subpopulations 

For a multipopulation model the choice of pa- 
rameters including predation rates can be enor- 
mously large. As a consequence more freedom is 
provided for choosing equilibrium densities that can 
fit the model ((l). Since it is impossible to inves- 
tigate the entire range of biologically plausible pa- 
rameters we make a particular choice of parameters 
that allow an easier comparison of the case n > 2 in 
(JTJ) with the model ([2]) . The details of the procedure 
are provided in Appendix. 

In this section we show the results of the nu- 
merical simulation for the model with n > 2 prey 
subpopulations. The system ^ for the case n > 2 
is integrated numerically. For the calculation of the 
stability diagrams at different fixed values of enrich- 
ment and transfer rate we perform 300 simulations. 
The results of the simulations for n — 2, 3, 4 and 7 
subpopulations are illustrated in Fig. [5] The per- 
centage of stable 3-species coexistence is calculated 
for every point in the parameter space with fixed 



enrichment K and limiting value t r . The colorcode 
is assigned according to the percentage of stable co- 
existence solutions found for 300 food webs. In all 
the replicas of the simulated system the steady state 
densities for G, F and Si are fixed (see (|A.9|) in the 
Appendix). Thus only the variations among pos- 
sible equilibrium densities {Sk}k^£i are examined. 
The constraints for the parameters of high dimen- 
sional system (TTJ are given in eqs. (IA.10|) - (IA.12[) in 
the Appendix. 

The stability diagrams in Fig. IA.5I show some 
similarities to the regions of coexistence in Figs. [5] 
and [7] found for the n = 2 subpopulation model 
The size of the stability zone expands with the 
increase of the transfer coefficient. At low transfer 
rates no stable persistence is found, but different al- 
ternative traits. The percentage of stable food webs 
with 3-species is substantially lower for a large sys- 
tem with n — 7 subpopulations than for n = 2,3. 
This reduction in stability is independent on the 
number of simulated food webs and a choice of main 
parameters of the system. It is possible that an in- 
crease of food web connectivity in this case impacts 
negatively the system (fT]) stability. Another feature 
is that for n = 7 the percentage of stable equilib- 
ria at a fixed enrichment value decreases for large 
values of t r unlike in previous cases in Fig. [5]a-c. 

The results of the numerical simulation demon- 
strate that for n — 2 subpopulation up to 95% 
of stable systems are found at a higher transfer 
rate and an intermediate enrichment. Second, for a 
larger food web with n — 7 subpopulations a higher 
percentage of stable steady states (up to 45%) are 
identified at low transfer rate and at high enrich- 
ment. 



We compare the results of simulation for four 
cases (n = 2, 3, 4, 7) at fixed enrichment K = 1.341 
and variable transfer coefficient in Fig. [9] The 
yields are derived for 1000 simulations of food webs. 
The estimations of the number of steady states 
show that for food webs with n < 4 a higher per- 
centage of solutions with a stable coexistence are 
identified than for food web with n = 7 pools. In- 
deed, the yield for n = 2 reaches almost 95% mean- 
while the percentage of stable food webs found for 
n = 7 saturates at 23% for large t r . The non- 
monotonic variations of the yields in Fig. [5] reveal 
a highly sensitive behaviour of the IGP model (fTJ) 
to a change in transfer rate in all cases. For n = 7 
the percentage of stable food webs reaches 41% at 
a low transfer rate. It decreases substantially for 
higher values of the transfer rate. For n < 4 there 
is an overall incline from 60% at t r ~ 0.1 to 95% at 
t r ~ 0.3 of stable configurations. 

Two types of stable equilibrium solutions are il- 
lustrated in Fig. [TUJ Both solutions are obtained 
inside the stable coexistence region as indicated in 
Fig. [5] The system (fTJ) is simulated with n = 5 
number of subpopulations and initial conditions as 
defined in the Appendix. For the steady state in 
Fig.[TU]a and the oscillatory state in Fig.[TU]b most 
of resource subpopulations are unstable and their 
densities rapidly decline to zero after some initial 
transient. Nevertheless, coexistence in the system 
is typically supported by one or two resource pools 
with non-zero densities. 
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Figure 6: Regions of stable (unstable) equilibria are marked 
according to species composition: 1(2) - stable (unstable) co- 
existence of resource and both predators G and F; 3 - exclu- 
sion of IG prey; 4 - exclusion of IG predator. The predation 
rates are: (a) / = 0.1, a = 0.0155 ;(b) / = 0.01, a = 0.0155; 
(c) / = 0.2, a = 0.0155; (d) / = 0.1, a = 0.065. The re- 
maining parameters are: r = 0.5,6 = 0.1, a' = 0.8, g = 
0.07,/' = 0.2, 9 ' = 0.5, m g = 0.02, m/ = 0.04 




Figure 7: Regions of stable (unstable) equilibria are marked 
according to stable trophic configurations: 1(2) - stable (un- 
stable) coexistence; 3 - exclusion of IG prey; (4) exclusion 
of IG predator. The predation rates are: (a) / = 0.1, a = 
0.0155 ;(b) / = 0.01, a = 0.0155; (c) / = 0.15, a = 0.065; 
(d) / = 0.1, a = 0.065. The remaining parameters are: 
r = 0.3, b = 0.02, a' = 0.8, g = 0.07, /' = 0.2, g 1 = 0.5, m g = 
0.02, = 0.04 
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Figure 8: The percentage of 3-species stable coexistence 
found for the model {TJ. The stability region is presented 
versus enrichment and transfer rate. Simulations for the 
four cases are given: (a) n = 2, (b) n = 3, (c) n = 4, 
(d) n = 7. Parameters are: r = 1, / = 0.1,6 = 0.1, a' = 



0.8, g = 0.07, /' = 0.2, g' = 0.5, m g 
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Figure 10: Two stable solutions for the system Q with 
n = 5 subpopulations. The equilibrium densities and in- 
teraction rates are described in Appendix. Parameters are: 
r = 0.7, / = 0.1,6 = 0.1, a' = 0.8, g = 0.07,/' = 0.2, g' = 
0.5, m g = 0.02, mj = 0.04, K = 1.8, t T = 0.2 



4. Discussion 
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Figure 9: The percentage of 3-species stable coexistence for 
the model N a i m = 1000 simulations are performed for 
the cases: n = 2, 3, 4, 7. Parameters are as in Fig. [8] Enrich- 
ment: K = 1.341. 



There is growing evidence from theoretical and 
empirical studies that creating additional trophic 
links have a stabilizing effect on food webs (jMoore , 



2005 



Ives and Carpenterl . 120071 ) . Generalized mod 



els reveal that the stability of food webs can 
be enhanced when species at higher trophic lev- 



els gr aze upon multiple prey species (jGross et al 



2009) 



In particular, for low dimensional food 
webs it is demonstrated that an addition of alter- 
native food resources can stabilize the interactions 



(jHolt and Huxell . 



20071) and open up a possibility 
for feedbacks on population dynamics due to ap- 
parent competition. The predictions of our model 
confirm the main conclusions give n in a theoretical 
study of an extended IGP model (|Holt and Huxell . 



20071 ). In the alternative formulations used here 



the IG prey has the access to an extra resource be- 
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yond the shared resource for which both predators 
compete. This extra resource is a more attractive 
resource item for the IG prey and is thus attacked 
at higher rates by the IG prey whereas the attack 
rate of the IG predator stays the same. Moreover 
the IG predator indirectly stimulated the growth of 
the IG prey population by providing this extra re- 
source. Our predictions tested by the application 
of a stability analysis are robust in the sense that 
they are independent of the form of the interac- 
tion term that is responsible for the availability of 
an additional resource. We demonstrate that for 
different formulations of the basic IGP model with 
the embedded interactions a stable 3-species coex- 
istence is ultimately reached whenever a moderate 
strength of the omnivorous links is used. 

Howe ver, the prob l em to relate the experimental 



findings (jLoder et al 



xxxxl ) to the theoretical pre- 



dictions of our model ^ still remains open. Since 
the experiment is aimed to observe a short term 
populations development it is not easy to find a di- 
rect correspondence between empirical population 
dynamics and theoretically predicted behaviour. In 
the 3 day batch culture experiment of Lodcr et al. 
with all 3 species present both predators G. dom- 
inans and F. ehrenbergii displayed positive growth 
while the prey population S. trochoidea displayed si- 
multaneously a sharp decline to almost zero. How 
can this behaviour be classified according to our 
theoretical model? Could it be a part of an oscil- 
latory cycle or an unstable state? It is not easy 
to answer these questions, however, we can make a 
guess that the short term evolution observed in the 
experiment recasts as a part of an oscillatory cycle 
for a periodic equilibrium state found at intermedi- 



ate immobilization. A similar type of experiments 
performed for various initial species densities could 
furnish a justification of this hypothesis. 

The above results demonstrate that a persis- 
tence of IG predator, IG prey and resource is 
achieved even at a low value of immobilization rate. 
Moreover, a significantly higher percentage of ob- 
served stable configurations is found when the im- 
mobilization and transfer links in ([1]) and @ are 
strengthened. Because our model is an oversim- 
plification of the experimental behaviour the parti- 
tioning of the parameter space according to stable 
versus unstable coexistence could be used as an ap- 
proximation of the population dynamics found in a 
real experimental situation. Firstly, the conditions 
for long term stable coexistence found by numerical 
simulations are not so easy to examine experimen- 
tally because of technical and te mporal restr i ctions . 



Experimental samples in ref. ([Loder et al 



are taken during 3 days of incubation due to a de- 
cline of the prey population. Secondly, due to the 
existence of stable limit cycles as predicted by our 
linear stability analysis (see Fig. ITU1) the oscillatory 
solutions go through a period of very low densities 
and might be driven to extinction in the presence 
of random fluctuations of the environment. 

We point out that our numerical simulations of 
the extended 3-species IGP module ([I]) do not ex- 
plore the entire swath of parameters and configura- 
tions for the steady states. Rather our analysis fo- 
cuses on explaining conceptual features of the IGP 
model with diverse prey populations. 

Earlier studies focussing on dynam- 
i cs of complex ecologi cal communities 
(jDeWitt and Langerhansl 120031 ) demonstrated the 
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importance of multiple prey traits in mitigating 
predator selection pressures and altering predator- 
induced behavioural shifts in natural environments 



(Lima and Dill 



1990: 



Trussell et al 



20021) . Our 



model can also be adapted to food web communi- 
ties in which differentiation among prey individuals, 
namely, variation in individual traits such as fitness 
and mobility, is a result of heterogeneities in their 
natural habitat and/or adaptation of the species to 
the local conditions of the habitat. We show that 
an existing diversity of resource items traits can sig- 
nificantly alter the emergent community patterns. 
Adding new subpopulations of resources with dis- 
tinct traits that are more vulnerable to an attack 
from the IG predator facilitates the coexistence 
of both IGP-related predators which compete for 
the common food resource. Thus, the presence of 
an alternative resource indirectly induces shifts in 
exploitative competition. 

It is important to note that a general math- 
ematical model with density- dependent interac- 
tions and immobilization do not render a unique 
theoret ical descript i on for the results of the exper- 



iment (Loder et al 



xxxxh . Our model predictions 



can be tested against alternative formulations. In- 
deed, the main features of the experimental system 
can be examined by the inclusion of predation rates 
that are dependent on the mobility of the resource 
species. Since slow and immobile individuals can 
also be found among mobile species one can use an 
inhomogeneous distribution of velocities of the re- 
source species in a theoretical model. To guarantee 
more benefit for an intermediate consumer in catch- 
ing a certain type of individual distinct predation 
rates should be assigned according to different ve- 



locities of resource species. Another question is if 
the growth rate of the IG predator will be affected 
by the inclusion of the time of resource capture. 
How will the inclusion of the time lag change the 
predictions of our immobilization model? These ex- 
tensions of a general IGP model will be a topic for 
our future investigations. 

Finally, we point out that it is of potential inter- 
est for biological control and conservation manage- 
ment to understand functioning of omnivory and 
IGP systems in relation to global changes of the 
environment. Since IGP food webs are widespread 
in natural communities their adaptation and re- 
silience behaviour is principal for understanding the 
restructuring of natural communities. 

5. Conclusions 

We have used three formulations of a general 
IGP model to explore the effects of increasing di- 
versity in the prey population on higher trophic lev- 
els. The reformulated IGP model alters the results 



from the basic IGP theory ( 


Polis anc 


Holt, 


1992 


Holt and Polis, 


1997; 


Diehl and Feifiel 


2000 


). We 



show that an increase of a number of trophic inter- 
actions in the system via differentiation of resource 
can stabilize the population dynamics of the IGP 
module. This conclusion holds for the densities of 
the IG prey that level up even when the IG preda- 
tor is a superior competitor for the common basal 
resource. 

The results of our numerical simulations can be 
summarized as follows. 

First, we show that for the system with the 
immobilization term up to three regions of stable 
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trophic configurations are observed along the en- 
richment gradient. Meanwhile at low enrichment 
both IG prey and IG predator are excluded, at high 
enrichment the presence of only small concentra- 
tion of immobilized cells is sufficient to facilitate 
the coexistence of the competitors in the IGP rela- 
tionship. Moreover the percentage of all admissible 
trophic configurations for the 3-species persistence 
inclines substantially for higher immobilization. 

Second, given that immobilization is high 
enough it prompts the exchange between pools of 
mobile and immobilized resource and facilitates fast 
decline of the mobile population and a growth of im- 
mobilized subpopulation. Meanwhile the exchange 
between the basal mobile resource and the preda- 
tors gets weaker due to the low density of mobile 
species the immobilized individuals become a ma- 
jor food resource for the predators. Because the IG 
prey is a superior competitor for immobilized re- 
source a robust coexistence of both predators will 
be easily supported. In addition, along an increas- 
ing gradient of immobilization the relative abun- 
dance of IG prey becomes higher than the abun- 
dance of IG predator. 

Restructuring of the basic IGP module by 
adding individual-to-individual turnover facilitates 
the coexistence and stabilizes the otherwise unsta- 
ble system. Moreover a strengthening of the inter- 
action link leads to a significantly broader range of 
enrichment values at which stable coexistence could 
be found. At low transfer rate two types of equilib- 
ria are observed: (i) if the IG predator is a superior 
competitor for the resources than at low enrichment 
both predators are excluded and at high enrichment 
only the IG predator stays in the system; (ii) an in- 



crease of the attack rates of the IG prey depresses 
the population of the IG predator until it is com- 
pletely excluded. 

Numerical simulations of food web ([TJ with 
n = 2,3,4 and 7 distinct pools demonstrate that 
the high dimensional food webs overall manifest far 
less stable behaviour than the food webs with only 
two distinct subpopulations. An interesting feature 
is that the percentage of stable states for n = 7 
substantially decreases from 40% to 23% with an 
apparent increase of the value of transfer rate. By 
contrast, for food webs with n < 4 an increase in 
transfer rate leads to the growth of the percentage 
of stable coexistence solutions from about 60% to 
95%. 

Appendix A. 

In the Appendix we review the steady state so- 
lutions for the Lotka-Volterra models ([2]) and ([3]) 
and provide Jacobian matrices to examine their lo- 
cal stability for the coexistence of both predators 
and the resource. 

First, equilibrium solutions are derived for the 
3-species model @ with zero immobilization (i m — 
0) and zero initial size of immobilized population 
(Si — 0). Second, the steady states are given for 
the model @ with immobilization (i m ^ 0). At 
last, the equilibrium solutions are presented for the 
system with the resource turnover ((3J). For every 
case various trophic configurations are considered: 
(i) exclusion of both predators, (ii) exclusion of IG 
prey or IG predator and (in) the 3-species coexis- 
tence. 
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Appendix A.l. Steady state solutions for 3-species 
model 

The equilibrium solution of @ for the 3-species 
coexistence without immobilization is stated as fol- 
lows: 



S eq = I /■ - 



amf rngf 



99 



K g \g> 
G eq = (m f - ff'S eq ){g'g)-\ 

F eq = (a'aS eq -m g )g~ 1 , (A.l) 

where g,g' > 0. The necessary condition for the 
coexistence requires that the right hand side is pos- 
itive in (|A.ip . 

The expressions for the equilibrium densities for 
the survival of the IG prey and the resource with 
exclusion of the IG predator at F eq = yield: 



^ = ^^ q = r -(l- S -k ). (A ,>, 



The condition for persistence of the IG prey and 
the resource reads: a'aK > m g . 

At zero density of the intermediate consumer 
(G eq = 0) one yields the steady states of the re- 
source S^ q and the IG predator F G : 

Sf,, 



(A.3) 



The densities are positive if and only if the condi- 
tion ff'K > TO/ holds true. 

Appendix A. 2. Steady state solutions and linear 
stability analysis 

Appendix A. 2.1. Model with immobilization 

Here we describe alternative steady state solu- 
tions and discuss their local stability derivation. 



Also the Jacobian matrix for the 3-species coex- 
istence is provided in the explicit form. 

For the model ([2]) with immobilization (i m ^ 0) 
we define a set of equilibrium densities to satisfy 
the equalities below: 

Qx = [r(l - SmK- 1 ) - aG 

-{f + im)F]S m , 
Q.2 = i m FS m -Si(bG + fF), 
Q 3 = (a'aS m + a'bSi - gF - m g )G, 



(A.4) 



Qi = (f'f(Si + S m )+ggfG-m f )F. 



The system (|A.5[) has four alternative solutions : (i) 
exclusion of both predators at S = K; (ii) exclu- 
sion of the IG predator (S^ ll G F ) at F = 0; (Hi) 
the coexistence of resource and the IG predator 
(S£;,Sf ,F G ) at G = 0; (iv) the 3-species coex- 
istence (S^,Sf,G e ,F e ). 

In the absence of F the immobilization mecha- 
nism is not active and the model ([2]) reduces to the 



syste m without immobilization (jDiehl and Feifiel 



EL 

iDic 



2000) where the equilibrium solutions written as 



(|A.2|> . Upon exclusion of the IG prey in (IA.5I) one 
obtains expression for the equilibrium densities of 
resource and the IG predator: 



/'(in 



/ + h 



f) 



1 



oG 

K 



'm Q G 



(A.5) 



Note that the size of mobile population S G is pro- 
portional to the size of immobilized population Sf* . 
As is expected the population of immobilized preys 
is impacted positively by the increase of immo- 
bilization. Although the predation pressures are 
equal for both resource subpopulations the immo- 
bile population extinguishes faster than the mo 
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bile population S^. Indeed, with the increase of 
predation rate / the following approximations hold: 
Sg ~ l/(i m + /) and Sf ~ l/f(i m + /). 

The equilibrium densities for the 3-species coex- 
istence are derived from (IA.5[) by setting the right 
hand side to zero. 

At last, to evaluate the stability of the equilib- 
rium solution S^, Sf, G e and F e one solves for the 
eigenvalues of the stability matrix : 



Di 

i F e 





D 2 
ba'G e 



-aS° 
-bSf 
D 3 



ffF e ff pe gg ,pe 



*m>Sm — fSl 

-gG e 



The matrix diagonal is written in terms of the 
equilibrium densities Sf^ , Sf , G e and F e as follows: 



D 



K 



S e m ,-bG e - fF e , 0,0 



(A.6) 



The solution (Sf n , Sf , G e , F e ) is globally 
asymptotically stable in the phase space 
(jSvirezhev and Loeofetl 119831 ) if the condition for 
stability is satisfied. For the stable coexistence it 
is necessary that the real parts of all four eigen- 
values Xi, (i — 1, . . .4) of the stability matrix are 
non-positive. To obtain the boundary for stability 
regions in the parameter space the eigenvalues of 
the above stability matrix are evaluated numeri- 
cally at different parameters combinations. The 
resulting stability diagrams are presented in Fig. [5] 
and Fig. [5] 

Appendix A. 3. Model with prey-to-prey interac- 
tions and n — 2 subpopulations 
As in the previous case the system [3] for n = 2 
pools permits four steady states: (i) the exclusion 



of the predators at Sx + S% = K; (ii) the exclusion 
of the IG predator (Sf , Sf , G F ) at F = 0; (m) the 
exclusion of the IG prey (Sf,Sf ,F G ) at G = 0; 
(iv) the coexistence of 3-species (Sf , Sj, G e , F e ). 

The solution for the coexistence of the resource 
and the IG prey in the absence of the IG predator 
is expressed as follows: 

ra'b 

F 



Sf = Kll- 



S. 



m g — a'aS{ 
aTb 



, G F = ^S[. (A.7) 



Note that at t r — the IG prey is excluded and 
steady state density for the resource approach the 
carrying capacity limit K . The positive solution of 
(|A.7p exists if the parameters satisfy the inequality: 



ra'b > t r m g . 



The steady state for the resource and the IG 
predator in the absence of IG prey yields: 

t r m 



r m f 

TO/ 



S 2 G = (I KI^^L-K. 
F G = — \ r (A.E 



Kt r ( r f m f 

rf V U f'f 



A nontrivial solution for the 3-species coexistence is 
found by solving for the equilibrium (Sf , Sf, G e , F e ) 
of the following system: 

= r(l - S^- 1 ) - aG e 

-fF e - t r S%, 
= t r Sf - bG e - fF e , 
= a'aSf + a'bS c 2 - gF e - m g , 
= /'/OS? + Sf) + g'gG e - m f . (A.9) 

Finally, the condition for the stable coexistence is 
provided by solving for the eigenvalues of the sta- 
bility matrix : 
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-r/KSf -t r Sf ~aSf -fSf 
t r S e 2 -bS e 2 -fS e 2 

aa'G e ba'G e -gG e 

ffF e ffF e gg ,pe q 

The eigenvalues of the stability matrix are func- 
tions of the constant rates (see Table. [2]) and the 
equilibrium densities Sf,S%,G e and F e . The con- 
dition for the stable coexistence requires that the 
real parts of the eigenvalues of the stability matrix 
are non-positive. 

Appendix A. 4- Model with prey-to-prey interac- 
tions and n > 2 subpopulations 

Here we define initial conditions that are used 
for the numerical simulation of the model (JTJ) with 
n > 2 number of subpopulations. Unlike for the 
system (J3]) with n = 2 subpopulations the equilib- 
rium densities can no longer be determined analyt- 
ically. To simplify the search for the equilibria in 
(TTJ) we implement several assumptions. 

We chose the parameters and equilibrium den- 
sities to fulfil several constraints provided below. 
Initial densities are equal to the steady states (|A.9|) 
for n = 2 subpopulations, namely: Si = Sf,S2 — 
5|, G = G e and F — F e . For the sake of simplic- 
ity the values for the remaining densities {Sk}^Z 2 
are defined to be less than the equilibrium density 
/Sf ■ Provided that the interaction rates {ck}k^i are 
randomly assigned values not exceeding t r /(n — 1) 
the equilibrium density S n can be found from the 
constraint: 

n 

E CjSj = t r Sl (A.10) 

We chose the zero rates of decline mi and mi for 
the first two subpopulations S\ and 5*2. The val- 



ues for the grazing rates {fk}k>3, are randomly 
assigned from the interval [0,f/n]. In the next 
step an antisymmetric matrix Q of individual-to- 
individual interactions is defined with the upper di- 
agonal coefficients {qjk}j<k that obey the inequal- 
ities: —0.1 t r < qj^k < 0.1 t r . The remaining lower 
diagonal coefficients {qkj}j<k should satisfy the an- 
tisymmetry relation: qj k — —qkj- The set of equa- 
tions in (fT]) for k = 3, . . .n holds true if the pos- 
itive mortality rates {mk}k>2 are expressed from 
the equations as follows: 

m k = CfcSi - E QjkSj - b k G - fkF. (A.ll) 

The equation (|A.11|) for k = 2 is used to solve 
for the attack rate f 2 . Finally, the attack rates 
WkVk=2 an d the feeding rates {f'k^Zl, are ran- 
domly assigned from the intervals [0,a'b/(n — 1)] 
and [0, f'f/(n — 1)] correspondingly. The choice 
of the feeding rates enables to reduce the preda- 
tion pressure on populations {Sk}^Z 2 by a factor 
of l/(n — 1). The remaining grazing rates /{ and 
b[ are defined from the relations: /{ = /'/ and 
b[ = a' a. Finally, two constraints hold to solve for 
the coefficients b' n and f' n : 

n n 

E f'i S i = f'f S ^ E h> S 3 = a ' bS 2- ( A - 12 ) 

3=2 j=2 

We want to emphasize that the attack rate b[ 
should be higher than any of the rates {&'-}™_ 2 . 
Nevertheless since a > b relation holds the total 
predation pressure of G summed over alternative 
pools exceeds the predation exclusively on S\ . This 
far the positive population density of the IG prey is 
maintained due to the consumption of alternative 
resources {Sj}j^ 2 - 



We provided a special assignment of parame- 
ters for the predation, feeding and mortality that 
enables to fulfil the condition of positive equilibria 
that are comparable to the realistic biodensities. 
Moreover the stability results can be conclusively 
compared with two different formulations ([T]) and 
©■ 
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